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We consider the cosmology of a "3-brane universe" in a five dimensional (bulk) 
space-time with a cosmological constant. We show that Einstein's equations 
r> I admit a first integral, analogous to the first Friedmann equation, which governs 

c^ • the evolution of the metric in the brane, whatever the time evolution of the 

metric along the fifth dimension. We thus obtain the cosmological evolution in 
the brane for any equation of state describing the matter in the brane, without 
needing the dependence of the metric on the fifth dimension. In the particular 
case p = wp, {w = constant), we give explicit expressions for the time evolution 
of the brane scale factor, which show that standard cosmological evolution can be 
obtained (after an early non conventional phase) in a scenario a la Randall and 
Sundrum, where a brane tension compensates the bulk cosmological constant. 
We also show that a tiny deviation from exact compensation leads to an effective 
cosmological constant at late time. Moreover, when the metric along the fifth 
dimension is static, we are able to extend the solution found on the brane to the 
whole spacetime. 
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1 Introduction 

There has been recently a lot of activity on the possibility that we live in a three- 
dimensional world embedded in a higher dimensional space. Modifying the old 
Kaluza-Klein picture 0, where the extra-dimensions must be sufficiently com- 
pact, these recent developpements are based on the idea that ordinary matter 
fields could be confined to a three-dimensional world, corresponding to our ap- 
parent Universe, while gravity could live in a higher dimensional space |^. The 
usual constraints on Kaluza-Klein models could therefore be relaxed 0], and 
"large" extra-dimensions would be conceivable, thus leading to a fundamental 
Planck mass much lower than its apparent three-dimensional value, even as low 
as the TeV scale. In the same spirit, Gogberashvili 0] and Randall and Sun- 
drum 0,0, reviving an older idea [0, have recently pointed out that the extra 
dimension (single in their case) need not even be compact. 

A first important question concerning these kinds of models is to check whether 
they give back standard four- dimensional gravity in the brane assumed to rep- 
resent our world. When the gravitational field generated by the brane itself is 
taken into account, it is not so easy to check that usual gravity is recovered. 
Randall and Sundrum showed that this is the case because of the existence of 
massless gravitons trapped in the brane, although small corrections are expected 
from massive modes due to the presence of a fifth dimension. Several groups have 
also begun to work on possible experimental signatures of these kinds of models 

i- 

In addition to the study of the weak field limit, one would expect that cos- 
mology as well could provide constraints on these types of models. An interesting 
aspect of cosmology in such models is that, thanks to the symmetries, it is pos- 
sible to deal with the non-linearities of the Einstein's equations and to detect 
effects that may not appear in a linearized treatment. The first crucial question 
is the viability of such models with respect to the cosmological evolution of our 
Universe. In a previous paper p, three of us have shown that, generically, the 
equations governing the cosmological evolution of the brane will be different from 
the analogous Friedmann equations of standard cosmology (the same holds for a 
brane with finite thickness |llO|)- Essentially, the difference lies in the fact that 
the energy density of the brane appears quadratically in the new Friedmann equa- 
tions in contrast with the linear behaviour of the usual equations. In the simplest 
model where the bulk is supposed to be empty (and stabilized), we found explicit 
solutions which showed that the corresponding cosmological scenario would be 
incompatible with nucleosynthesis constraints. 

Following our work, it was noticed by two groups ([|r^,|[12|) that a possible way 



to reconcile brane cosmology with the required standard cosmological scenario (at 
least since nucleosynthesis) would be to use Randall and Sundrum's approach 0, 
i.e. to put a (negative) cosmological constant in the bulk and to introduce in the 
brane, in addition to ordinary matter, a constant tension that would exactly 



compensate the bulk cosmological constant so that the, quadratic, leading order 
due to the brane contribution would be cancelled, leaving the next, linear, order. 
These works however only provided an approximate perturbative analysis of 
the cosmological evolution. Now, exact solutions of the global Einstein's equa- 
tions are extremely useful, and until now, such solutions with a bulk cosmological 
constant exist only when the spacetime is static (0), or in the extensions given 
by Kaloper [^ (see also IQ and |[l^), where the brane tension does not ex- 



actly compensate the bulk cosmological constant, so that dynamical evolution is 
possible, but in the simplest form of exponential expansion. Some other exact 



solutions are also known in the context of general supergravity |16| or Horava- 



Witten supergravity |T^. General conditions for static universes can be found in 



IS 



The purpose of the present work is to solve the five-dimensional Einstein's 
equations for any type of matter in the brane with a cosmological constant in the 
bulk. We first show in section 2 that Einstein's equations admit a first integral, 
which in particular provides directly the cosmological evolution of the brane. In 
the following section, we show that an additional assumption, namely that the 
metric along the fifth dimension does not evolve in time, enables us to solve for 
the whole space-time metric, i.e. to find explicitly the dependence of the metric 
on the transverse coordinate as well as time. Finally, in the last section, we 
solve analytically, for different cases, the new Friedmann equation obtained in 
the present work and discuss the consequences. 

2 Solving Einstein's equations 

Let us present here the general framework. We shall consider five-dimensional 
spacetime metrics of the form 

ds'^ = QABdx dx^ = g^ydx^dx'^ + b'^dy'^, (1) 

where y is the coordinate of the fifth dimension. Throughout this article, we will 
focus our attention on the hypersurface defined hj y = 0, which we identify with 
the world volume of the brane that forms our universe. Since we are interested 
in cosmological solutions, we take a metric of the form 

ds'^ = —n^ij, y)dT'^ + a^ij, y)'~^ijdx^dx^ + b^ij, y)dy'^, (2) 

where 'jij is a maximally symmetric 3-dimensional metric {k = —1,0,1 will 
parametrize the spatial curvature). 

The five- dimensional Einstein equations take the usual form 

Gab = Rab — -^Rqab = f^ Tab, (3) 



where Rab is the five-dimensional Ricci tensor, R = q^^Rab the scalar curvature 
and the constant k is related to the five-dimensional Newton's constant G(5) and 
the five-dimensional reduced Planck mass M(5), by the relations 

k' = 87rG(5) = M^l (4) 

Inserting the ansatz (^ for the metric, the non-vanishing components of the 
Einstein tensor Gab are found to be 
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In the above expressions, a prime stands for a derivative with respect to y, and 
a dot for a derivative with respect to r. 

The stress-energy-momentum tensor can be decomposed into two parts, 

-^ B^T Blbulk + -^ Slbranc; (9) 

where T^ ^Ibuik ^^ ^^^ energy momentum tensor of the bulk matter, which will be 
assumed in the present work to be of the form 

T^ sLuik = diag (-PB, Pb, Pb, Pb, Pt) , (10) 

where the energy density ps and pressures Pb and Pt are independent of the 
coordinate y. Later, we will be specially interested in the case of a cosmological 
constant for which —pB = Pb = Pt- 

The second term T^ BLianc corresponds to the matter content in the brane 
{y = 0). Since we consider here only strictly homogeneous and isotropic geome- 
tries inside the brane, the latter can be expressed quite generally in the form 

^^Blb.anc = -^dmg{-pb,Pb,Pb,Pb,0) , (11) 

where the energy density pb and pressure pb are independent of the position inside 
the brane, i.e. are functions only of time. We may also add some other brane 
sources with similar energy momentum tensor. 



The assumption that T05 = 0, which physically means that there is no flow 
of matter along the fifth dimension, implies that 6*05 vanishes. It then turns 
out, remarkably, that the components (0, 0) and (5, 5) of Einstein's equations (|^) 
(with d^) and (§)) in the bulk, can be rewritten in the simple form 

(12) 
(13) 
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ka^. (14) 

Since T° q = —pB is here independent of j/, one can integrate ([12|), which 
gives 

F+—a''pB + C = Q, (15) 

o 

where C is a constant of integration which a priori depends on time. Assuming 

in addition that T° q = T^ 5, one finds using the time derivative of (0) and 

the y— derivative of (pTSj) that pB is constant in time. This also implies that C is 

constant in time. In order to deal with the last component of Einstein's equations 

(see (0)), it is convenient to use the Bianchi identity. 



V^G^o = 0, (16) 

which can be rewritten (using Gqs = 0, and T^ q = T^ 5) in the form 

dr (^) = laa'l^Gi. (17) 

One finds (using (jT^)) that this equation is identically satisfied if —pB = Pb- 
Hence, when the bulk source is a cosmological constant, any set of functions a, 
n, and h satisfying (|15D or, more explicitly. 
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together with Gqs = 0, will be solution of all Einstein's equations (^, locally in 
the bulk. 

The brane can then be taken into account by using the junction conditions 



IS], which simply relate the jumps of the derivative of the metric across the brane 



to the stress-energy tensor inside the brane. This procedure is described in detail 



in our previous paper [^ in the context of five-dimensional brane cosmology (see 
also Ijl^)- The relevant expressions are 

'"'' - 4p. (19) 
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where the subscript for a, b, n means that these functions are taken in y = 0, 
and [Q] = Q{0^) — Q{0^) denotes the jump of the function Q across y = 0. 

Assuming the symmetry y ^^ —y for simplicity, the junction condition (p^) 
can be used to compute a' on the two sides of the brane, and by continuity when 
y — s> 0, ([l8|) will yield the generalized (first) Friedmann equation (after setting 
no = 1 by a suitable change of time): 

The salient features of this equation are that, first, the bulk energy density enters 
linearly, second, the brane energy density enters quadrat ically, and finally the 
cosmological evolution depends on a free parameter C (related to the choice of 
initial conditions in the whole space-time), whose infiuence corresponds to an 
effective radiation term (from the standard point of view, i.e. interpreting linearly 
the additional term)Q. One can note that this equation is a generalization of the 
exact equation obtained in our previous paper 0. This equation is enough to 
study the cosmological evolution in the brane, independently of the metric outside 
and in particular of the time evolution of b. The analysis of this equation will be 
postponed until section 4. In the next section we will be interested in obtaining 
an explicit solution for the whole five dimensional metric. 

3 Explicit dependence on the fifth dimension 
for a stabihzed bulk 

In this section, it will be shown that, with the help of an additional assumption, 
namely that the fifth dimension is static, in the sense that 

6 = 0, (22) 



^From a technical point of view, the presence of a free parameter, in contrast to the standard 
cosmological equations, can be explained in the following way. In four-dimensional cosmology, 
the component (0,0) of Einstein's equations yields a first integral, because it contains only first 
derivatives in time. In the present work, we have used the component (0,0) to integrate in y, 
whereas the time integration is operated with the help of the component (5,5) which contains 
second derivatives in time. 



it is then possible to solve the full space-time metric, i.e. to determine the explicit 
dependence of the metric on the coordinate y. The restriction (^) allows us to 
go to the gauge 

h=l. (23) 

It then follows immediately from equation G05 = 0, that n can be expressed in 
terms of a according to the relation 

- = a(t), (24) 

n 

where a is a function that depends only on time (and not on y). Inserting this 
into (|l^) yields the following differential equation: 

a^ + k- {aa')' = —psa^, (25) 

o 

which is valid everywhere in the bulk (but not in the brane) on the two sides of 
the brane separately. It can be integrated in y, yielding 

a^ = A cosh(/xy) + B sinh(/iy) + C, (26) 

with 

/ 0TI2 

(27) 

(28) 
(29) 



(30) 

for pb = 0. In the following, we will focus on the first case ps < 0, but all the 
equations will apply as well to the case ps > 0, up to the transformation p -^ ip, 
B -^ -iB. 

The coefficients A, B, C, D, E are functions of time, and C is expressible in 
terms of a as 

C = 3^, (31) 

and the others can be determined by the junction conditions. The symmetry 
y ^^ —y imposes the relations A^ = A_ = A, B^ = —B_ = B between the 
coefficients on the two sides of the brane. Using ([T9|) and (0), one then finds 

pB K^ 2pB 2 / 1 \ ,^^, 
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Note that one can check exphcitly energy conservation in the brane from these 
relations, i.e. 



ao 



0. 



(33) 



Let us finally fix our temporal gauge, by imposing the condition uq = 1, i.e. 
we choose as time what corresponds to the cosmic time in the brane. The function 
a is then simply do (see equation (|2^) ), and (|^) specialized to y = gives 
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with, from equations (^TJ) and ([21|), 
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2 V 6pB, 
Using the first relation in (p2D, the remaining coefficient is determined as 

B = -'^al (36) 



3/i " 
Substituting back the coefficients A and B thus obtained in the general expression 



(^61) and using (|2l|) , one finally obtains the following expression for the scale factor 
a{t,y): 
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The other coefficient of the metric, n{t,y), follows directly from the above ex- 
pression with the help of the relation (El), i.e. 



n{t,y) 
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(3^ 



Therefore, from the two functions ao(^) and pb{t) restricted to the brane, which 
are obtainable by solving the coupled system constituted of the two equations 
(I2TD and (^3|), one is able to infer the extension of the metric in the bulk away 
from the brane. 

Let us mention briefiy the case pb = 0. Following the same procedure as 
above, one would end up with the following expression for the scale factor: 
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(39) 



In the special case C = 0, one recognizes the hnear solutions obtained in our 
previous paper 0. 

We have not specified here whether the fifth dimension is compact or not. If 
it is compact, then one must check that one gets in the end a globally well defined 
solution (see 0). 

The next section will be devoted to the determination of the evolution of the 
scale factor on the brane, aQ(t), using equations (|2ll) and (|33|), which - we wish 
to emphasize - does not depend on the hypothesis (p2D of the present section. 



4 Cosmological scenarios 



In this section we will examine the consequences of equation (21) on the cos- 
mological evolution inside the brane. In general, one still finds the quadratic 
behaviour in the brane energy density which leads to a non conventional cosmol- 
ogy 1^. However, as suggested by Jill ^^^ BI1I> *^ "^i^^ show that one can indeed 



implement the idea of Randall and Sundrum ^ in the cosmological context to 
recover standard cosmology. We will find exact solutions that match the early 
non conventional cosmology to the subsequent ordinary cosmology. 

Let us thus assume that the energy density in the brane can be decomposed 
into two parts, 

Pb = pA + P, (40) 

where p^ is a constant that represents an intrinsic tension of the brane and p 



stands for the ordinary energy density in cosmology. Substituting in (^TD one 
gets 

If we follow Randall and Sundrum by choosing pA such that 

yPB + '^PK = 0, (42) 

then one sees that standard cosmology is recovered with the identification ( [ITTI] , 
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SvrG ^ ^, (43) 



when p ^ Pa- Let us however keep pA unspecified at this stage. In the case p = 
and C = one would recover the solutions given by Kaloper [p^ . 

To get analytic solutions we will now assume that ordinary matter is described 
by an equation of state of the formp = wp, with w constant. Then, using equation 



(33) for p and p, one can write 

p = p*(ao/a*)~^ g = 3(l + u;), (44) 



where p* and a* are constants. Inserting this expression in pl|), it is possible to 
integrate exphcitly the resulting equation in the case where C = and k = 0. 
Assuming the first term on the right hand side of ( Jl) to be positive and defining 
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one finds, for A > 0, 
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(the origin of time being chosen so that ao(0) = 0). It is clear, from the latter 
expression, in the case A = 0, how one passes from a very early universe, charac- 
terized by a non-conventional evolution a{t) ~ t^^"^, to a late time phase described 
by standard cosmology, a{t) ~ t^/'^. In the case A > 0, with A sufficiently small, 
one obtains three successive phases, a non-conventional phase dominated by p^, 
a conventional phase dominated by p and, finally, an exponential phase, where 
A plays the role of an effective cosmological constant in our Universe. Therefore, 
brane cosmology is able to produce a late time acceleration phase such as seems 
to be needed from the latest cosmological observations (pO[), at the price how- 
ever of a tiny mismatch in the compensation between p^ and the tension of the 
brane. The required fine-tuning expresses the brane cosmological version of the 
well-known cosmological constant problem. 

Let us finally examine the case of a universe filled with radiation {w = 1/3). 
It is then possible, still with A; = 0, to integrate explicitly the case C ^ 0, because 



the C-term has the same dependence on qq as the term proportional to p in (^l]) . 
Moreover, the free parameter C can be constrained by nucleosynthesis. Indeed, 
at the time of nucleosynthesis, the universe is dominated by the radiation energy 
density, which can be written 



TT 



PraditN) — g*l^T^, 



(48) 



where g^, is the effective number of relativistic degrees of freedom at that time. 
In the standard model, g,,{standard) = 10.75, and any deviation A^f,, is strongly 
constrained by the observed abundances of hght elements, typically Agf^, < 2. 
In our model, since the additional C-term evolves like radiation, it can also be 
seen effectively as additional relativistic degrees of freedom, subject to the usual 
constraint, so that 
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at nucleosynthesis. 

To conclude, we wish to recall the main results of this work. First, it has 
been shown that one can obtain a first integral of Einstein's equations, which 
provides, on the brane, a relation analogous to the (first) Friedmann equation 
and which depends only on the geometry and matter content of the brane, except 
for a constant parameter. Second, when 6 = 0, one can extend explicitly the 
solution found on the brane to the whole spacetime. Finally, we have shown, ex- 
hibiting exact solutions, that brane cosmology appears compatible, at late times, 
with standard cosmology, in a Randall-Sundrum type model. In this respect, an 
important question for future investigations is wether an (exact or quasi) Randall- 
Sundrum configuration can be naturally reached in a dynamical scenario for the 
very early universe. 



Note added: Several papers which discuss related issues have appeared quasi 



simultaneously to this one [21, B2] (see also 
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